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Abstract. The freeze out (FO) of the expanding systems, created in relativistic heavy ion collisions, is 
discussed. We start with kinetic FO model, which realizes complete physical FO in a layer of given thickness, 
and then combine our gradual FO equations with Bjorken type system expansion into a unified model. 
We shall see that the basic FO features, pointed out in the earlier works, are not smeared out by the 
expansion. 



PACS. 51.10.+y - 24.10.Nz - 25.75.-q 

I At highest energies available nowadays in relativistic 
. heavy ion collisions at RHIC the total number of the pro- 

■ duced particles exceeds 6000, therefore one can expect 
, that the produced system behaves as a "matter" and gen- 
' erates collective effects. Indeed strong collective flow pat- 
I terns have been measured at RHIC, which suggests that 

■ the hydrodynamical models are well justified during the 

■ intermediate stages of the reaction: from the time when lo- 
, cal equilibrium is reached until the freeze out (FO), when 
' the hydrodynamical description breaks down. During this 
I FO stage, the matter becomes so dilute and cold that par- 
I tides stop interacting and stream towards the detectors 

■ freely, their momentum distribution freezes out. The FO 
, stage is essentially the last part of a collision process and 

the main source for observables. 

Nowadays, FO is usually simulated in two extreme 
ways: A) FO on a hypersurface with zero thickness, B) 
FO described by volume emission model or hadron cas- 
cade, which in principle requires an infinite time and space 
for a complete FO. At first glance it seems that one can 
avoid troubles with FO modeling using hydro+cascade 
two module model [1], since in hadron cascades gradual 
FO is realized automatically. However, in a such a scenario 
there is an uncertain point, actually uncertain hypersur- 
face, where one switchs from hydrodynamical to kinetic 
modeling. First of all it is not clear how to determine such 
a hypersurface. This hypersurface in general may have 
both time-like and space-like parts. Mathematically this 
problem is very similar to hydro to FO phase transition on 
the infinitely narrow FO hypersurface, therefore for exam- 
ple all the problems discussed for FO on the hypersurface 
with space-like normal vectors will take place here. An- 
other complication is that while for the post FO domain 
we have mixture of non-interacting ideal gases, now for 
the hadron cascade we should generate distributions for 



the interacting hadronic gas of all possible species, as a 
starting point for the further cascade evolution. The vol- 
ume emission models are based on the kinetic equations 
[5113] defining the evolution of the distribution functions, 
and therefore these also require to generate initial distri- 
bution functions for the interacting hadronic species on 
some hypersurface. 

In this work we present FO model which allows us to 
study FO in a layer of any thichness, i, from to oo, 
and which connects the pre FO hydrodynamical quanti- 
ties, like energy density, e, baryon density, n, with post 
FO distribution function in a relatively simple way. Many 
building blocks of the model are Lorentz invariant and can 
be applied to both time-like and space-like FO layers. In 
this work we are going to include Bjorken like expansion 
in our FO model, in contrast to the older versions P',T,'51[S1 
U]. In this latter case the FO layer is a domain restricted 
by two hypersurfaces r = ti and t = ti + L (t is the 
proper time). 

We are going to review briefly all the steps done to 
derive our model. We will skip all the detailed derivations, 
refering to the corresponding publication, and will show 
and discuss only a small part of results, due to limited 
space. 

Starting from the Boltzmann Transport Equation, in- 
troducing two components of the distribution function, 
/: the interacting, /% and the frozen out, ones, (/ = 
P + f^)j and assuming that FO is a directed process (i.e. 
neglecting the gradients of the distribution functions in 
the directions perpendicular to the FO direction compar- 
ing to that in the FO direction) we can obtain the follow- 
ing system of the equations [51[5] : 

ds Tpo Tth ' ds Tpo 
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The FO direction is defined by the unit vector dcr^. FO 
happens in a layer of given thickness L with two paral- 
lel boundary hypersurfaces perpendicular to da^, and s — 
dafj,x'^ is a variable in the FO direction. We work in the 
reference frame of the front, where dcr^j is either (1, 0, 0, 0) 
for the time-like FO, or (0, 1,0,0) for the space-like FO. 
The tfo is some characteristic length scale, like mean free 
path or mean collision time for time- like FO. The rether- 
malization of the interacting component is taken into ac- 
count via the relaxation time approximation, where fi 
approaches the equilibrated Jiittner distribution, feq{s), 
with a relaxation length, Tth ■ The system ([T]) can be solved 
semi-analytically in the fast rethermalization limit [6] . 

According to the above references, the basis of the 
model, i.e. the invariant escape probability within the FO 
layer of the thickness L, for both time-like and space-like 
normal vectors is given as [6l[7ll9] 



L 



where is a particle four-momentum, is the flow ve- 
locity. In fact the model based on the escape rate ([2]) is a 
generalization of a simple kinetic models studied in Refs. 

which can be restored in the L ^ oo limit. Here 
we will concentrate on the time-like case only, where the 
above function is unity. 

Simple semianalytically solvable FO models studied in 
[3ll4l[5l[6l[7] are missing an important ingredient - the ex- 
pansion of the freezing out system. The open question is 
whether the features of the FO, found in those papers, will 
survive if the system expasion is included. In this work 
we present a model which includes both gradual FO and 
Bjorken-like expansion of the system. And we will see that 
the answer is "yes" - the basic features of the post FO dis- 
trubitions will not be smeared out by the expansion. 

First, let us remind the reader the basics of the famous 
Bjorken model. Bjorken model is one-dimensional in the 
same sense as discussed before eq. JT]) - only the proper 
time, T = yW^^x^, gradients are considered. Here the 
reference frame of the front, da^ = (1, 0, 0, 0), is the same 
as the local rest frame, u'^ = (1,0,0,0). The evolution 
of the energy density and baryon density is given by the 
following equations: 



de e + P dn n 
dr T ' dr t 



(3) 



where P is the pressure. The initial conditions are given 
at some t — tq. 

Applying our FO model to such a system, we obtain: 



df{r') = -^-L^Pir') + ^ [Ur') - f{r')] 

TFO L-r Tth 



d//(r')- + 



dT' L 



(4) 
(5) 



Tpo L-t' 

where FO begins at r = ri and r' = r — ri . Taking the fast 
rethermalization limit, similarly to what is done in [5], we 
can obtain simplified equations for which is a thermal 



distribution /'(r) — feq{T), as well as for e*,n* and 



de' 
d? 

def 
d^ 



e* L dn' 
Tfo L — t' ' dT' 

e' L dn^ 



TFO L~ t' 

n' L 



(6) 
(7) 

Tfo L — t' ^ dT' tfo L — t' 

Now the idea is to create a system of equations which 
would describe a fireball which simultaneously expands 
and freezes out. Let us put our two components (e = 
e* A- e^) into the first equation of ([3]) and do some sim- 
ple algebra: 



L 



e-* e 

r Tfo L ~ t' tfo L ~ t' ^ 



de' def 
dT dT 

(8) 

where last two terms add up to zero; the free component, 
of course, has no pressure. So far our eq. ([8]) is completely 
identical to the first equation of ([3]). Our assumption is 
that our system evolves in such a way that eq. ([8]) is satis- 
fied as a system of two separate equations for interacting 
and free components [10] : 



de' 
'd? 



P' 



L 



def 



TFO L + Tl-T 
3' L 



(9) 



(10) 



T TFO L + Tl-T 

Similarly we can obtain equations for baryon density |10j : 



dn' 
dnf 



n 



n 



L 



T Tfo L -\- Ti — t 
n^ L 



(11) 



(12) 



T Tfo L + Ti — t 

Thus, finally, we have the following simple model of 
fireball created in relativistic heavy ion collision. 
Initial state, r = tq: e(To) = cq, n{To) = no. 
Phase I, Pure Bjorken hydrodynamics, tq < t < ti 



e('r) = eo (-^^ , n{T) = no (J^^ 



(13) 



where P = c^e - equation of state (EoS) in general form. 
Phase II, Bjorken expansion and gradual FO, 

Ti < T < Ti + L 



eo 



(?)' 



'TO 

no\- 



L + Tl-T 



L + Tl 

L 



L/tfo 



L/tfo 



(14) 



(15) 



With these last equations we have completely determined 
evolution of the interacting component [10] . Knowing e' (r) 
and EoS we can find temperature, Ti{T). Due to symmetry 
of the system uf(T) = u^(to) = (1,0,0,0). Finally, /'{t) 
is a thermal distribution with given Ti{T), n'{T), u'^{t). 
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However for us the more interesting is free component, 
which is the source of the observables. Eqs. ()9ll0p give us 
the evolution of the e / and n f , and one can easily check 
that these two equations are equivalent to the following 
equation on the distribution function: 



= 5.0 fm/c, X, =0.5 fm/c 
1 fo 



dT 



L 



TFO L + Tl-T 



(16) 



The measured post FO spectra are given by f^{L + ti). 

Aiming for a qualitative illustration of the FO process 
we show below the results for the massless ideal gas with- 
out conserved charges with Jiittner equilibrated distribu- 
tion (P' = eV3, = ^T^. r(T,|p|) = ^e-|Pl/^'M). 
We have taken the following values of the parameters: 
To = 0.05 /m, r,(To) = 835 MeV- n = 5 /m, T,(ti) = 
Tpo = 180 MeV\ tfo = 0.5 fm and we present results 
for different values of FO time L. 
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Fig. 1. Evolution of the temperature of the interacting mat- 
ter for different FO layers. Ti(ro = 0.05 fm) = 835 MeV, 
Tfo ~ 180 AleV. "No Freeze Out" means that we used stan- 
dart Bjorken hydrodynamics even in phase II. 



Fig. [T] shows the evolution of the temperature of the 
interacting matter. As it was already shown in iSjiTj the 
final post FO particle distributions are non-equilibrated 
distributions, which deviate from thermal ones particu- 
larly in the low momentum region. By introducing and 
varying the thickness of the FO layer, L, we are strongly 
affecting the evolution of the interacting component, see 
Fig. [11 but we gain see the universality of the final post FO 
distribution: for L > 2tfo it already looks very close to 
that for an infinitely long FO calculations: see Fig. O The 
inclusion of the expansion into our consideration does not 
smear out this very important feature of FO. More results 
can be found in [TU] . 

In our opinion these results may justify the use of FO 
hypersurface in hydrodynamical models for heavy ion col- 
lisions, but with a proper non-thermal post FO distribu- 
tions. If the FO layer is thick enough, say L > 2tfo, then 
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Fig. 2. Final post FO distribution for different FO layers as 
a function of the momentum in the FO direction, p = in 
our case {p^ = = 0). The initial conditions are specified 
in the text. "No Expansion" curve is given by the analytical 
expression [5]: f^{p) = ~j^Ei{-~^). 



it doesn't matter how thick was FO layer, we do not need 
to model the FO dynamics in details. Once we have a good 
parameterization of the post FO spectrum (still asymmet- 
ric, non-thermal), for example analytical post FO distri- 
bution obtained in Ref. [5] (see caption of Fig. [2]), then 
the parameters of this distribution can be found from the 
conservation laws, as it is usually done for sharp FO, with 
some volume scaling factor to effectively account for the 
expansion during FO. It is important to always check the 
non-decreasing entropy condition p!Ol[TT] to see whether 
such a process is physically possible. 
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